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ABSTRACT 

We investigate the diffusive siiock acceieration in the presence of the non-resonant streaming 
instability introduced by Bell (2004). The numerical MHD simulations of the magnetic field 
amplification combined with the analytical treatment of cosmic ray acceleration permit us to 
calculate the maximum energy of particles accelerated by high- velocity supernova shocks. The 
estimates for Cas A, Kepler, SN1006, and Tycho historical supernova remnants are given. We 
also found that the amplified magnetic field is preferentially oriented perpendicular to the shock 
front downstream of the fast shock. This explains the origin of the radial magnetic fields observed 
in young supernova remnants. 

Subject headings: cosmic rays- acceleration of particles-sliock waves-supernova remnants- instabilities 



1. Introduction 

The instabilities produced by energetic parti- 
cles are important phenomena accom panyin g the 
diffusive sh ock ac celerat ion (K rymsky (|l977l ): Ax- 
ford et al. (Il977l) : Bell (|l978l ). Blandford and Os- 
triker ( IQTsh ) of cosmic rays in supernova rem- 
nants (SNRs). The scattering of energetic par- 
ticles both upstream and downstream of a su- 
pernova shock is supplied by magnetic inhomo- 
geneities existing in the shock vicinity. It was 
suggested that this may be the result of a res- 
onant streaming instability that develops due to 
the presence of diffus ive streaming of accelerated 
particles (Bell (119781 )). The total magnetic field 
may be amplified if the energy of unstable magne- 
tohydrodynamic (MHD) waves becomes compara- 
ble with the energy of the background magnetic 
field. 



Such an amplification seems quite possible be- 
cause quasilinear theory of the resonant streaming 
i nsta bility allows it (see e.g. McKenzie & Volk 
(|M)). 

The presence of amplified magnetic fields in 
young SNRs is well established now. It was indi- 
cated by the radio observations of SNRs and was 
attributed to the Ray leigh- Taylor instability of the 
contact discontinuity between the gas of supernova 
ejecta compressed at the reverse shock and the cir- 
cumst ellar g as compressed at the forward shock 
(Gull (|l975l )'). The radial magnetic fields inferred 
from measurements of the radio-polarization in 
young shell type SNRs (see e.g. Milne (1987.) ) may 
appear as this instability develops. 

However, the discovery of thin X-ray filaments 
coinciding with the position of the forward shock 
in young gal actic S NRs (Gotthelf et al. (l200lh . 



Hwang et al. (|2002h . Vink & Laming (|2003h . Long 
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etal (|2003f) . Bamba et al. (|2003l) . Bamba et al. 
([20051)) has led to the conclusion that the mag- 
netic field is amplified just at the forward shock. 
This conclusion does not depend on the nature 
of the mechanism which produces such filaments: 
the fast synchrotron cooling of electrons acceler- 
ated at the forwa rd sho ck (see e.g. the analysis by 
Berezhko et al. (|2002l )). or the dissipation of the 
MHD turbulence and the corresponding decrease 
of t he ma gnetic field, as was suggested by Pohl et 
al. (|2005h . 



Recently Bell (|2004l ). using the dispersion re- 
lation for c oUisio nless MHD waves derived by 
Achterberg |l983'), found a new regime of a non- 
resonant streaming instability. In the presence of 
the strong electric current of accelerated particles 
a non-oscillatory purely growing MHD mode ap- 
pears at spatial scales smal ler than the gyroradius 
of the particles. Bell (|2004l ) performed MHD sim- 
ulations and showed that the magnetic field may 
be strongly amplified. 

This instability is investigated in more d etail i n 
our companion paper (Zirakashvili et al. (|2008l ). 
Paper I) via numerical MHD simulations. Here 
we combine these simulations with the analyti- 
cal treatment of diffusive acceleration at a plane 
steady-state shock. It allows us to estimate the 
maximum energy of the accelerated particles and 
to obtain the value of the amplified magnetic field. 

The paper is organized as follows. The ana- 
lytical model of diffusive acceleration at the plane 
shock is considered in Sect. 2. The MHD simula- 
tions are described in Sect. 3 and 4. The maximum 
energies of accelerated particles are estimated in 
Sect. 5. Sect. 6 contains the discussion of obtained 
results. The summary is given in the last Sect. 7. 

2. Acceleration at the plane parallel shock 

We shall consider the generation of MHD tur- 
bulence and the particle acceleration in a simple 
one-dimensional case and assume a steady state in 
the reference frame of the shock. The applications 
of our results to real three-dimensional shocks are 
considered in the next Sections. 

The upstream plasma moves with a velocity u = 
Ml from —oo along the z axis. The plasma velocity 
downstream u — U2 — ui/a drops by a factor of 
a at the shock front located at z = 0. Here a is 
the shock compression ratio. We shall consider a 



parallel shock; therefore the mean magnetic field 
Bq is in z direction. 

We shall also neglect the effects of the mean 
electric field Eg directed along the mean magnetic 
field. The electric field modifies the cosmic ray 
transport equation (see Paper I). The isotropic 
part of the cosmic ray momentum distribution 
N(j), z) obeys the following cosmic ray transport 
equation upstream and downstream of the shock: 



d ^ dN ^dN _ ^ 
dz " dz dz 



(1) 



Here is the parallel diffusion coefficient of the 
energetic particles. The cosmic ray distribution 
N{p) is normalized as = / 47rp^iV(p), where 
the cosmic ray number density. 
The function N is continuous at the shock. The 
boundary condition of the cosmic ray flux conser- 
vation at the shock front, z — can be written 
as 



7s dp dz 



- D, 



ON 
dz 



(2) 



where 7^ = 3ct/((T — 1), No{p) is the distribution 
function at the shock, Di and D2 are the parallel 
diffusion coefficients upstream and downstream of 
the shock, respectively. 

We impose an additional boundary condition 
N — Oa,tz = —L. This qualitatively describes 
the escape of the highest energy particles from a 
SNR with the distance L being of the order of the 
supernova shock radius R. 

The solution of Eq. (1) in the upstream region 
z < may be written as 



1 — exp J uidzi/ Di{zi,p) 

N{z,p)^No{p) ^ . (3) 

1 — exp J uidzi/ Di{zi,p) 

In the downstream region z > the solution is 
simply N = Nq. The boundary condition (2) gives 
the ordinary differential equation for Nq{p): 



dp 



IsNoip) 



1 — exp — J uidz/Di{z,p) 



(4) 



2 




0.01 



Fig. 1. — The momentum distribution No{p) at 
the shock (solid Hne) and the cosmic ray flux J{p) 
(dashed line) at the absorbing boundary aX z = 
—L. The shock compression ratio is cr = 4. 



where the function b{z) is given by the expression 



b{z) 



tt/2 



(3/8) sin^ Ode 



J d^kByyCk, z)S{kz cos 9 + smt 



(6) 

Here Byy (k) is the spectrum of the y-component of 
the random magnetic field. For the isotropic ran- 
dom field this function is 6 = / d'^kBisotr{k)/ {2k). 
The solution of Eq.(4) can then be written as 
No{p) (X p~^^nQ{p/pm) where the function no(s) 
with the argument s = p/pm describes the shape 
of the spectrum in the cut-off region: 



no(s) = exp 



dsi/si 



exp(s]^ 



(7) 



Since the non-resonant instability produces a 
random magnetic field with scales smaller than 
the gyroradius of the particles, the smal l-scal e ap- 
proximation of Dolginov and Toptygin (!l967 ) can 
be used for the calculation of the scattering fre- 
quency v which determines the diffusion coeffi- 
cient, see Paper L The corresponding mean free 
path A = v/iy is proportional to the square of the 
particle momentum. 

If particles are scattered by the small-scale 
isotropic field, the scattering frequency does 
not depend on the pitch-angle 6 and the diffu- 
sion coefficient along the mean magnetic field 
is Z?|| = v^/3v. The scattering frequency v — 

j0 J d3fcBi,otr(fc)/fc (cf. Paper I) is deter- 
mined by the spectrum of the isotropic random 
magnetic field Bisotiik). It is normalized as 
(SB'^) = / Bisou {k)d^k, where (SB"^) is the mean 
square of the random magnetic field. 

The scattering frequency depends on the pitch- 
angle in a more general case when the random 
field is statistically isotropic in the plane that is 
perpendicular to the mean field direction (see Ap- 
pendix A). Let us introduce the momentum pm 
defined as: 



2 37r q Ui 

Prn o 



dzh{z), 



(5) 



It is convenient to write down the distribution 
Nq{p) in terms of the cosmic ray energy flux Fe 
at the absorbing boundary at z — —L. This is the 
energy flux of the highest-energy particles escap- 
ing from a SNR (the so-called run-away particles). 
It may contain an essential part rjesc = ^Fe/ pu\ of 
the kinetic energy flux in particular when 

the acceleration is efficient and the shock structure 
is modified by the pressure of accelerated particles. 
Here p is the plasma density. 

The momentum distribution at the shock front 
can then be written as 



iVo(p) 



VescPUi 



4— 7s 



-no{p/p„ 



(8) 



8ttcp'^'P„ 

Here / = dss^-'^^nois){eyip s^^ - 1)-^. The 
function No{p) and the flux of run-away particles 
at the absorbing boundary J{p) — ui7Vo(p)/(exp(p^/p^) 
1) are shown in Fig.l. 

We use such a normalization of the spectrum of 
accelerated particles mainly because the parame- 
ter rjesc is directly related to the number density 
of accelerated particles in the cut-off region. The 
electric current of these particles drives the non- 
resonant instability (see below). It is possible to 
use other parameters instead of rjesc, e.g. an in- 
jection efficiency of the thermal ions at the shock 
front, but since its relation with the high-energy 
end of the spectrum is not so straightforward, we 
prefer to use rjesc- 
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In addition we have a physical reason to use this 
normahzation. Since generally the shock propa- 
gates in the medium with a very large diffusion co- 
efficient (e.g. it is higher than 10^^ cm^s~^ in the 



interstellar medium, cf. Berezinskii et al. ([I99O0), 
the accelerated particles are effectively scattered 
only near the shock front where the level of the 
self-excited MHD turbulence is high. The accel- 
erated particles with the maximum energies run 
away from this region to the outer space. Thus the 
acceleration at the plain shock with the absorbing 
boundary considered in this paper simulates the 
real situation when the acceleration by the three- 
dimensional supernova shock is considered. 

The possible effect of the shock modification by 
the cosmic ray pressure was not taken into account 
in Eq.(4). The modification may strongly change 
the spectrum of particles and make it c oncav e at 
low energies (see e.g. Malkov & Drury (j2001r ) for 
a review). However, even in this case the cut- 
off of the spectrum is described accurately by Eq. 
(7). The total shock transition that includes the 
precursor created by the pressure of accelerated 
particles and the thermal sub-shock may be con- 
sidered as a sharp discontinuity for the cosmic ray 
particles with the highest energies. 

The calculation of the diffusive electric current 
of accelerated particles with the use of Eq. (8) 
gives 



jd{z) 



2cpml 



exp 



exp(s ^) — 1 



no(s), 



(9) 



where a(z) — J dzib{zi). 

-L 

We shall use the values 7s = 4 and 7 = 1/4 
below. This corresponds to an unmodified strong 
shock. However, the accepted value of 7^ does not 
strongly differ from 7^ ~ 3.7 — 3.8, typical for the 
moderately modified shocks. 

3. Modeling of the non-resonant instabil- 
ity with diffusive shock acceleration 

We can now model the magnetic field amplifica- 
tion in the vicinity of the shock which accelerates 
particles. We shall seek the steady state solution 
for the spectrum of accelerated particles and the 
time-averaged MHD spectra. 
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Fig. 2. — The dependence of the normalized max- 
imum momentum on the normalized shock ve- 
locity vl^ . The results were obtained using the res- 
olution 64'^ (crosses), 128"^ (open circles) and 256'^ 
(stars) . The analytical approximation and the up- 
per limit according to Eq. (19) are shown by the 
solid and dashed lines, respectively. 



Since the shock velocity is much higher than 
the phase velocities of MHD waves and the tur- 
bulent velocities of the plasma upstream of the 
shock, we can model the dependence of the MHD 
spectra on the distance from the shock front via 
the investigation of the temporal evolution of the 
instability in the simulation box moving with the 
gas flow in the direction of the shock front. This 
strongly reduces the size of the simulation box in 
the z-direction and permits to obtain the numeri- 
cal results with a good numerical resolution. The 
computation time is also significantly reduced. 

We shall assume that the shock propagates in a 
medium with the density po, the gas pressure Pq 
and the Alfven velocity Va = Bq/ y/AirpQ. 

The details of the numerical method were given 
in Paper I. 

The dimensionless time i, the space coordinate 
z and the velocity u are determined ds, t = tVakQ, 
z — koz, u — u/Va- Here is the wavenumber 
that corresponds to the real size 27r/fco of the box. 
The dimensionless density p and the electric cur- 
rent J can be expressed via the magnetic field -Bo 
and the Alfven velocity Va &s p — AirpV^ /Bq and 
J = iirj/ckoBo- 

The real size of the simulation box is small 
in comparison with the characteristic scale of the 
spatial distribution of the accelerated particles up- 
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stream of the shock. At i = the box is placed 
at z = —L, where the initial background ran- 
dom magnetic field corresponding to the isotrop- 

ically distributed Alfven waves with the one di- 
mensional spectrum cx and the amplitude 

1/2 

= 0.09i?o is preset. We use the gas adi- 
abatic index 7 = 5/8 and the parameter (3=1. 
Here (3 = 47rPo/Bo- 

The box moves with the mean flow speed ui 
towards the shock. The MHD equations which in- 
clude the Lorentz force produced by the electric 
current of accelerated particles are solved numer- 
ically in the three dimensions (see Paper I for de- 
tail). 

Let us fix the momentum Pm- At every instant 
of time the diffusive electric current of accelerated 
particles, that drives the instability, is calculated 
according to Eq. (9). It may be rewritten in di- 
mensionless units as 



00 

J{i) = JoJ 



exp (5s 2 di'bii') 

72 



exp(s 2) — 1 



Ms). (10) 



Here g is some arbitrary dimensionless constant. 
We change the integration over dzi in Eq. (9) to 
the integration over time dt' = dzi/ui in Eq. (10). 
The dimensionless current Jo can then be written 
in terms of the physical parameters as 



Jo 



2?7e 



1/2 

' n 



Sttc 



1/2 



(11) 



The simulation is performed up to the point in 
the dimensionless time im at which the value of the 
integral Jq (iF6(?) reaches g~^. This corresponds 
to the box arrival to the position of the shock. 
It means that we have found the size L for the 
given maximum momentum pm- Or, inversely, the 
momentum p^ can be found from Eq. (5): 
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Fig. 3. — Dependence of the physical quantities 
on z. Top: random magnetic field (solid line), 
turbulent velocity (dashed line) and sonic veloc- 
ity (dotted line). Bottom: momentum distribu- 
tion N{pm) (sohd line), diffusive electric current 
J (dashed line), scattering frequency v{pm) nor- 
malized to the gyrofrequency ^{pm) — <iBo/pm 
(dotted line) and the mean electric field Eq (thick 
solid line). The calculations were performed for 
the normalized shock speed u'^^ — 4900 km s~^. 

to the value about 4 and therefore the character- 
istic scale of magnetic field is smaller than the size 
of the box. 

It is convenient to present the numerical results 
in terms of the normalized shock velocity u\ de- 
fined as 



^ = iz} ( 



mc 



V 25c 



1/2 



qBoL 



(12) 



Several runs were performed to scan the broad 
range of physical parameters. We used the value 
Jo = 64 and several values of 5 > 1 and Jo = 
64(7~^ for g < I. This choice limits the character- 
istic wavenumber of the generated magnetic field 



Va 



U-t = Ul r 

^ VlOkms"^ 



-3/4 1 /„ 

V0.05/ 



(13) 



and the normalized maximum momentum p'^, 
given by 
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Fig. 4. — The one-dimensional spectra of the per- 
pendicular component of the random magnetic 
field B±{kz) (solid line), turbulent velocity U±{kz) 
(dashed line) and magnetic hehcity H{kz) (dotted 
line) at the shock front. The calculations were per- 
formed for the normalized shock speed u'-^^ = 4900 
km s~^. All spectra are normalized to the mean 
square of the perpendicular component of the ran- 
dom magnetic field. 



Pm =Pn 



Va 



10 km s" 



-1/2 



L 
1 pc 



i—Y 

(14) 

We normalize the shock velocity and maximum 
momentum using the parameter value ijesc = 0.05. 
This value corresponds in particular to the case 
when the energetic spectrum of CR particles at 
the shock front is a power-law with the expo- 
nent 7s — 2 = 2 and the total CR pressure equals 
0.5^0^*1 • The dependence of the normalized max- 
imum momentum p'^ on the shock velocity u[ is 
shown in Fig. 2. 

The spatial dependence of the random magnetic 
field, the turbulent velocity, the sonic velocity, the 
momentum distribution N{pm), the diffusive elec- 
tric current, the scattering frequency i'{pm) and 
the mean electric field Eq obtained for the nor- 
malized speed u[ = 4900 km s~^ are shown in 
Fig.3. 

The spectra of the magnetic field, the plasma 
velocity and the magnetic helicity at the shock 
front are shown in Fig. 4. 

For physical applications it is useful to know the 
probability distributions functions (PDFs) of dif- 
ferent physical quantities. The PDF of the random 
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Fig. 5.— PDF of the magnetic field Pb{B) (thin 
solid line), turbulent velocity Pu{u) (thin dashed 
line) and density Pp{p) (thin dotted line) at the 
shock front. The numerical calculations were per- 
formed for the normalized shock speed u'^ = 4900 
km s~^. 



magnetic field, the turbulent velocity and the gas 
density obtained in our simulations at the shock 
front are shown in Fig. 5. 

As it is seen in Fig. 3 the diffusive electric cur- 
rent sharply increases near the shock. As a result 
the perpendicular components of the turbulent ve- 
locity also increase and the kinetic energy of the 
random motions just upstream of the shock is an 
order of magnitude larger than the magnetic en- 
ergy. 

It is important that the accelerated particles 
are concentrated in the shock vicinity. This jus- 
tifies the use of the planar geometry even for real 
three-dimensional shocks. The amplification of 
the MHD turbulence takes place in the broader 
region where the diffusive electric current of run- 
away particles is not small. 

On the whole our modeling corresponds to the 
following physical picture. Let us consider a vol- 
ume element at some distance from the supernova. 
Shortly after the explosion the run-away particles 
reach the volume and drive the streaming instabil- 
ity. Our numerical modeling shows how the mag- 
netic fluctuations are amplified in the volume ele- 
ment. For simplicity we assumed that the elec- 
tric current was constant at all times after the 
explosion before the shock arrival. This assump- 
tion is strictly valid for a steady state plane shock 
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Table 1: Results of the modeling of the non-resonant instability with diffusive shock acceleration 



u[, 10^ km/s^ 


1.55 


2.76 


4.90 


8.25 


13.9 


23.3 


39.2 


p'^, TeV/cb 


9.46 


22.7 


56.4 


120 


227 


443 


854 


(SB^f' /Bo" 


1.36 


4.68 


12.6 


24.2 


44.8 


87.3 


187 


0.77 


0.71 


0.77 


0.85 


0.85 


0.84 


0.82 


{Su^f'/Va^ 
Cs/Vj 


5.15 


23.6 


49.9 


94.1 


166 


361 


808 


1.42 


3.70 


11.9 


31.6 


68.4 


139 


289 


(K/10 km s-i)V4(^^^jo.05)-i/2:/(p^)/r!(p^)g 
102(K/10 km s-i)-i^o/So'^ 
102(K/10 km s-^)-^'''q^/p,nd 
WD[pm)/uiV' 


1.44 


2.63 


3.95 


4.96 


7.21 


12.6 


22.4 


0.016 


0.14 


0.74 


2.86 


10.6 


43.7 


209 


0.348 


1.03 


2.03 


3.30 


5.00 


7.52 


11.0 


9.2 


6.8 


6.3 


6.3 


4.9 


3.3 


2.1 


{Va/IQ km s-i)-i/*(r?e.c/0.05)i/2A:oP™c/gBo' 


0.194 


1.09 


6.15 


14.6 


34.8 


82.7 


197 



"Normalized shock velocity ti'j^. 

''Normalized momentum p'^ of accelerated protons. 

'^Ratio of rms of the random magnetic field to the mean magnetic field. 

''Ratio of rms of the random parallel component of the magnetic field to the corresponding perpendicular component 

{{8Bl+5Bl)/2f' 

Ratio of rms of the turbulent velocity to the Alfven velocity. 
■'^ Ratio of the sound velocity to the Alfven velocity. 

fRatio of the scattering frequency y{pm) to the gyrofrequency Q{pm) = qBo/pm- 
''Ratio of the mean electric field to the mean magnetic field. 

'Ratio of the potential energy of the particle q(p corresponding to the electric potential to the energy PmC- 
''Ratio of the diffusion length Dijim) /u\ and L 

'Product of the wavenumber fco and the gyroradius pmc/qBo- The wavenumber determines the real size 27r/A,'o of the simulation 
box. 

Note. — All quantities from the lines c-k are calculated at the shock front. 



(see Fig. 3) and only qualitatively valid for three- 
dimensional shocks. The dimensionless time 
in our calculations corresponds to the supernova 
remnant age T. If the shock radius R increases as 
R cx t"^° , where mo is the expansion parameter, 
then T = moR/ui and we should use the relation 
L = moR = uiT for the parameter L in Eqs (12) 
and (14). Note, that the electric current increases 
according to our plane shock modeling when the 
shock comes close to the volume element. 

A summary of the numerical results obtained 
for different normalized shock velocities u'l is given 
in Table 1. 

4. MHD modeling in the shock transition 
region and downstream of the shock 

The density fluctuations arising during the de- 
velopment of the non-resonant instability, when 
the expanding magnetic spirals collide each other 



(cf. Paper I), play a crucial role in the shock tran- 
sition region. It is well known, that the propaga- 
tion of the shock in the medium even with small 
density fluctuations is accompanied by shock front 
distortions and the appearance of vortex plasma 
motions downstream of the shock (see e.g. Kon - 
torovic h (1959,) . McKenzie & Westphal (|l968h . 
Bykov (|l982l )). These motions can amplify the 
magnetic field downstream of the shock even in the 
absence of accelerated particles. This effect was 



recently observed by Giacalone & Jokipii (|2007f ) 
in their 2D MHD numerical simulations. The re- 
sults of low resolution three dimensional calcula- 
tions for the MHD evolution of an adiabatic su- 
pernova remnant in a non-uniform and turbulent 
inte rstella r medium were presented by Balsara et 
al. (|200lh . 

We investigated this phenomenon and studied 
the evolution of MHD turbulence in the down- 
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MAGNETIC FIELD 




Fig. 6. — Slices of the magnetic field strength 
(top) and gas pressure (bottom) through the cen- 
ter of the box obtained at i = 0.12 in the 
shock transition region and downstream of the 
shock. The scaling of the magnetic field strength 
is logarithmic between 12.6i?o (black) and 126i?o 
(white). The scaling of the gas pressure is also log- 
arithmic between IO^Bq/Att (black) and IO^Bq/Att 
(white). The calculations were performed for a 
shock propagating with a speed ui = 3000 km s~^ 
in the medium with Va = 10 km s~^. 

stream region. The simulation box with the size 
27r X 27r X Aw, stretched in the z direction was used. 

The gas flows with the shock velocity ui along 
the z-axis into the box from its left boundary. At 
the initial moment of time the flat shock front 
is located at z = 27r. The plasma compression 
and heating downstream are taken in accordance 
with the Rankine-Hugoniot conditions. During 
the simulation the plasma magnetic field, density 
and pressure are prescribed at the left boundary 
in accordance with the spatially periodic numer- 
ical solution, found in the previous Section. The 
plasma leaves the system at the right boundary, 
where the homogeneous boundary conditions are 
prescribed. The periodic boimdary conditions in 
the perpendicular directions are used. We per- 




Fig. 7. Dependence of the physical quantities 
on z in the shock transition region. Top: ran- 
dom z- and perpendicular magnetic field compo- 
nents (normal solid and thin solid lines respec- 
tively), z- and perpendicular components of the 
turbulent velocity (thick and thin dashed lines re- 
spectively). Bottom: sonic velocity (dotted line), 
mean gas density (dashed line), mean gas velocity 
(solid line) in z direction . The calculations were 
performed for the shock propagating with a speed 
Ui = 3000 km s~^ in the medium with Va = 10 
km s~^. 
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Fig. 8. — The one-dimensional spectra of the per- 
pendicular component of the random magnetic 
field B±{kz) (solid line), turbulent velocity U±{kz) 
(dashed line) and magnetic helicity H{kz) (dot- 
ted line) downstream of the shock. The calcula- 
tions were performed for a shock propagating with 
a speed ui = 3000 km in the medium with 
Va = W km s~^. All spectra are normalized to 
the mean square of the perpendicular component 
of the random magnetic field. 
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Fig. 9.— PDF of the magnetic field Pb{B) (thin 
solid line), turbulent velocity Pu(u) (thin dashed 
line) and density Pp{p) (thin dotted line) down- 
stream of the shock. The analytical approxima- 
tion for PDF of the magnetic field (23) (thick dot- 
ted line) is also shown. The calculations were 
performed for a shock propagating with a speed 
Ui = 3000 km s~^ in the medium with 14 = 10 
km s~^. 



formed 3D MHD simulation with 256^ x 512 cells 
using the numerical method described in Sect. 3. 
The electric current of accelerated particles was 
switched off. 

The obtained numerical results at the shock ve- 
locity ui = 3000 km s~^ and the Alfven velocity 
Va = 10 km s~^ are shown in Fig. 6 and Fig. 7. We 
used the numerical solution with the normalized 
shock velocity u'l = 4900 km s~^ described in the 
previous Section. The corresponding parameter 
rjesc = 0.14. 

The real size of the box in z direction = 
4:TT/ko is related to the distance to the absorbing 
boundary L as = 4nLVa/{uiim)- The value 
i„i = 1.02 corresponds to the numerical solution 
with the normalized shock velocity u'^ = 4900 km 
s~^ described in the previous Section. 

The slices of the magnetic field strength and 
the gas pressure in the YZ plane are shown in 
the top and bottom panels of Fig. 6 respectively. 
The strong distortions of the shock front and the 
sonic waves propagating in the downstream re- 
gion are clearly seen in the bottom panel. The 
shock front is shifted to the right due to the in- 
teraction with the fluid elements which have the 



enhanced density. The magnetic structures down- 
stream are stretched along the direction of the flow 
(top panel). 

The dependence of the plasma parameters 
averaged in the perpendicular XY plane on z- 
coordinate is shown in Fig. 7. Strong fluctuations 
of the plasma motions with amplitude of the or- 
der of one third of the shock velocity ui exist 
downstream. It is remarkable that these random 
motions occur mainly in z direction. The mag- 
netic field is also stretched in this direction. The 
z-component of the random magnetic field is a 
factor of 1.4 larger than the perpendicular com- 
ponents downstream of the shock. 

The spectra of the magnetic field, the plasma 
velocity and the magnetic helicity downstream of 
the shock are shown in Fig. 8. Note, that the mag- 
netic helicity is not zero in this region. 

PDF distributions of the random magnetic 
field, turbulent velocity and plasma density ob- 
tained downstream of the shock are shown in 
Fig.9. 

The numerically calculated magnetic compres- 
sion ratio (Tb is about 3.0 in this simulation with 
the shock compression ratio a close to 4. This 
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value of Us is close to one that is expected if the 
isotropic random magnetic field is compressed at 
the shock front: as = \/ (Scr-^ + l)/3. Since in 
the real situation the pressure of accelerated cos- 
mic ray particles should be taken into account, the 
shock compression will be higher. We found that 
the use of the adiabatic index 7 = 4/3 with the 
corresponding shock compression ratio about 6 re- 
sults in a magnetic compression factor as close 
to 3.5, that is smaller than the expected value 
Ob ~ 5. Probably this is due to the numerical 
dissipation of the field. It is necessary to have a 
better numerical resolution in the last case since 
the stronger compression of the field at the shock 
front results in the stronger numerical dissipation. 
We obtained the close value of ub = 2.6 in the 
simulation with 7 = 5/3 and the lower numerical 
resolution 256 x 128^. This demonstrates that the 
resolution 512 x 256^ is good enough for the sim- 
ulation of the shock with the compression ratio 4. 
We found that the corresponding values of ctb dif- 
fer significantly from each other in the case of the 
higher shock compression ratio. 

5. Analytical estimate of the maximum 
energy of accelerated particles in SNR 

Now we can estimate analytically the maximum 
particle energy achieved in the process of acceler- 
ation. We shall use Eq. (5) and assume that the 
random field is isotropic and concentrated at the 

wavenumber fc, so that b = B'^/{2k) where the rms 

1/2 

of the random magnetic field Br — {SB^') . Let 
us change the integration on dz to the integration 
on time dt = dz/ui in Eq. (5). The evolution of 
the amplified magnetic field can be approximately 
described as 



Br{t) 



Bb exp(7ma 

2So(l+7n 



,t), t < tl 

ax(i-ii)), t > tl 



^-1 ln2B!i 

/max Bb 



(15) 

Here Bt, is the initial value of the amplified field. 
At the first stage t <ti the field is amplified expo- 
nentially in time with the maximum growth rate 
7max = jdBo/2cpoVa (cf. Paper I). The corre- 
sponding wavenumber k is k — 2'Kjd/cBf) at this 
stage. At the second stage t > ti the field is am- 
plified linearly in time. The wavenumber k may 
be found from the expression kBr = A'Kjd/c. This 
dependence of random magnetic field on time is 
in a qualitative agreement with the simulations of 



the non- resonant instability (cf. Paper I). 
The integration of Eq. (5) gives 



37r q'uiB'^Va 



fmax^ 



l + (l+7max(im-tl))^, 



(16) 

Here tm = L/ui. The electric current jd that 
drives the instability is constant over the whole 
region — L < z < except the narrow zone near 
the shock (see Fig. 3). The input of this zone was 
neglected in Eq. (16). This assumption introduces 
a relatively small error into the calculation of Pm, 
since Eq. (16) is an implicit equation for pm- The 
growth rate 7max is inversely proportional to Pm 
(seeEq. (17) below). Thus remains in Eq. (16) 
only through 7inax after the parenthesis. This pro- 
duces an error in the evaluation of and p,,„ that 
is not larger than 10 percent. 

Using Eq. (9) we estimate jd = rjescqpul/ {2cpm) 
and 

r]esc qufBo 



7n 



4 VaC'^Pn 



(17) 



From these two equations we finally obtain the 
value of the amplified field 



Br = 2Bri 



1/4 (18) 
1*1 > u* 



and the maximum energy of accelerated particles 



In" 



PmC 



r]escqulBoL 



2Bo 



,Ui < u» 



Ml > U* 



2^ - 1 



1/41 



(19) 

where the velocity = {2AttcV^ /-q^g^y^^ . 

Using the first line of expression (19) for an 
arbitrary shock velocity ui we obtain the upper 
limit of the maximum energy. The growth rate 
(17) calculated with this maximum energy is high 
enough to provide the considerable magnetic field 
growth during the age of a sup ernova remnant (to 
be compared with Bell ( 20041 )'). The upper limit 
and the maximum energy (19) for In (2Bq/ Bh) = 5 
are shown in Fig. 2. We found that for ui » 
Eq. (18) gives the amplified field a factor of 2 
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smaller than the numerical results (see Table 1). 
It is because we neglected the growth of the dif- 
fusive current jd in the narrow zone just adjacent 
to the shock in our analytical estimates. However, 
this does not influence the estimate of the maxi- 
mum energy (19). 

For the fast shocks ui >> the amplified up- 
stream magnetic field may be estimated using the 
formula 



47r 



- . 2 Va 

— ^ 2rjescPUi\ — . 



(20) 



For velocities ui > u» Eq. (19) may be rewrit- 
ten as 



Prn^ 



21 TeV 



0.05 \^ 



1000 km s 



H ^pc 



In ( ^ ) - 1 



- 1 



1/4' 



(21) 

where Z is the charge number of accelerated parti- 
cles, U3 = 1730 km s~^, riH is the hydrogen num- 
ber density in cm^^ and Lpc is the size L in par- 
sees. Note that the denominator of this equation 
contains the normalized shock velocity (12). 

For three-dimensional shocks this equation may 
be considered as a rough estimate of the maximum 
energy of accelerated particles. Then the size L is 
related to the shock radius R and the remnant age 
T as L ^ rrioR = uiT . 

6. Discussion 

The initiation of the MHD streaming instability 
by accelerated particles in the shock precursor is 
an integral part of the efficient diffusive shock ac- 
celeration process. The instability is non-resonant 
if the electric current of accelerated particles is 
large enough, so that the normalized shock veloc- 
ity exceeds some critical value 



Ui > (4cKV'7e.c)'/' 



1340 km s" 



Va 



10 km s 



2/3 



/ r]esc \ 

V0.05/ 



-1/3 



(22) 

This equation is equivalent to Eq. (18) in Paper 
I. 



The non-resonant character of this instability 
means that the principal scale of the growing ran- 
dom magnetic field remains smaller than the gy- 
roradius of particles with momentum close to its 
maximum value pm (t) . The treatment of the non- 
magnetized particle scattering is relatively simple 
and it allowed us to fulfill the numerical simula- 
tion of the diffusive shock acceleration accompa- 
nied by the strong streaming instability and to 
make the corresponding analytical estimates. The 
calculated amplification of magnetic field in young 
SNR proved to be very significant. The maximum 
energy of accelerated particles is mainly limited by 
the finite time of the growth of magnetic distur- 
bances. The analytical estimate of the maximum 
particle momentum is given by Eq. (19) and it is 
in agreement with our numerical results presented 
in Fig. 2. The resulting maximum momentum is 
not as high as one may expect using the Bohm 
diffusion coefficient in the amplified field since the 
scattering by the small scale field is not so efficient. 

The small-scale field approximation is broken 
for particles with relatively low energies when 
their gyroradii in the amplified magnetic field 
are smaller than the principal scale of the field. 
Roughly it occurs when the scattering frequency 
of the particle is smaller than its gyrofrequency 
in the amplified field. The ratio of these quanti- 
ties for particles with the maximum momentum 
Pm is the ratio of the numbers from the 7th and 
3rd lines of Table 1 . This ratio is close to 1 for the 
smallest considered normalized velocity u[ = 1550 
km s~^ when the small-scale approximation is only 
marginally valid. The ratio decreases if the shock 
velocity increases. This ratio is close to 1/3 for 
the normalized velocity 4900 km s~^ that is the 
representative value for the young historical su- 
pernovae (see below) . This means that for this 
shock velocity the small-scale approximation is 
valid for particles with the normalized momentum 
larger than p'^/i. For particles with smaller en- 
ergies the amplified magnetic field can be consid- 
ered as the mean large-scale field. These particles 
are resonantly scattered by the magnetic inhomo- 
geneities from the inertial range of the magnetic 
spectrum (see Fig. 4) or by the magnetic perturba- 
tions produced by the resonant streaming instabil- 
ity of these particles (cf. Pelletier et al. ( 2006f )). 

It is interesting that formula (21) without the 
root in the denominator may be used to estimate 
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Table 2: Maximum energies and amplified magnetic fields in historical SNRs. 









nH" 


Bo^ 








PmC* 




PmC* 








Bd 












0.01 


0.05 


0.14 


0.01 


0.05 


0.14 


0.01 


0.05 


0.14 




Tycho 


435 


4500 


0.3 


5.0 


1360 


3040 


5090 


19 


76 


170 


21 


110 


260 


300 


SN1006 


1000 


4300 


0.1 


5.0 


860 


1920 


3220 


22 


100 


240 


8.4 


43 


120 


140 


Kepler 


400 


5300 


0.35 


5.0 


1670 


3800 


6360 


30 


115 


250 


33 


160 


350 


215 


Cas A 


330 


5200 


3.0 


10.0 


2220 


4960 


8300 


63 


230 


500 


120 


510 


980 


485 



"Age of SNR, yr 

'Forward shock velocity, km 

'^Hydrogen number density of the circumstellar medium, cm~^ 
''Magnetic field strength in the circumstellar medium, /iG 
■^Normalized shock velocity, km s~^ 
■'^Calculated maximum energy pmC, TeV 

s Calculated downstream magnetic random field strength, fiG 

''Downstream magnetic field determined from the thickness of X-ray filaments, 

Note. — The nu mbers from the columns b, c are taken from Vink ||2006| ). the numbers in the last column h are according to 
Volk et al. l|2005f ) 



the maximum energy of particles, accelerated at 
slow astrophysical shocks when the condition (22) 
is violated and the resonant streaming instabil- 
ity should be taken into account. The matter is 
that the expression for the increment of the non- 
resonant instability (17) differs from the increment 
of the resonant instability only by a factor of the 
order of unity (see e.g. Berezinskii et al. (1990)). 
This is why the expression (21) may be also used 
for slow shocks provided that the wave damping 
does not quench the development of instability 
that is relatively slow in this case. It may be in 
a highly ionized medium where the damping of 
MHD waves on neutrals is negligible. Nonlinear 
damping of Alfven waves may be also depressed 
in a plasma with (3 > 0.3, since the presence of 
waves, propagating in the opposite to the cos- 
mic ray streaming dire ction is necessary for this 
damping (Zirakashvili ('2000')). Such waves may 
appear only in a coUisionless low-(5 plasma due to 
the no nlinear induced scattering (Livshits & Tsy- 
tovich ( 1970f )). or due to the three- wave interac- 
tions in the framework of magnetohydrodynamics 
(Chin & Wentzel (|l972l) l. 

Because the particles at the end of the spectrum 
are scattered by the small-scale magnetic field, the 
spectrum has the universal shape in the cut-off re- 
gion described by Eqs (7) and (8) in the case of 



high velocity shocks. This is important for the cal- 
culation of gamma-ray production by the nucleon 
component in SNRs. 

The comparison of the predicted values of am- 
plified fields given in Table 1 with the values de- 
rived from the thickness of X-ray filam ents of 
young SNRs (see e.g. Volk et al. (|2005h ) shows 
reasonable agreement, if the acceleration efficiency 
is not low: rjesc > 0.05. 

The calculated downstream magnetic fields and 
the maximum energies obtained for the histor- 
ical SNRs with known ages are given in Table 
2. We use the same supernova parameters as 
Vink ( 2006f ). We assumed the standard interstel- 
lar value of the magnetic field Bq — 5 /iG for SNRs 
Kepler, Tycho and SN1006. The accepted mag- 
netic field strength Bq = 10 fiG for Cas A super- 
nova gives the value of the Alfven velocity about 
10 km s~^ that is a reasonable number for the stel- 
lar wind produced at the Red Supergiant stage 
of the likely progenitor of this supernova. The 
magnetic field amplified upstream of the shock is 
determined by the value of the normalized shock 
velocity u[ and is given in the 3rd line of Table 
1. This almost isotropic random magnetic field is 
further amplified by a magnetic compression fac- 
tor as = 4 in the shock transition region. 

One of the key parameters in our calculations 



12 



rjesc is in principle determined by the injection 
efficiency of thermal particles in the process of 
acceleration and by the degree of shock modifi- 
cation by th e cosm ic ray pressure, s ee Be rezhko 
and Ellison (|l999l) and Blasi et al. (|2005l ). The 
value of rjesc can not be calculated theoretically 
yet and we used different rjesc in our estimates. 
The results obtained for three values rjesc — 0.01, 
Vesc = 0.05 and rjesc = 0.14 are shown in Table 
2. The value rjesc = 0.14 is realized for the plain 
cosmic ray modified shock with the compression 
ratio (7 = 6 and the thermal sub-shock compres- 
sion ratio CTs — 2.5. The second value rjesc = 0.05 
corresponds to the situation when due to some rea- 
son the shock modification is weaker. The value 
Vesc = 0.01 corresponds to the non-modified shock 
with the cosmic ray pressure Pgr about ten per- 
cents of the ram pressure poUi- 

The maximum energies of particles accelerated 
in Kepler, Tycho and SN1006 supernovae that all 
are of the la type are about 100 300 Z TeV. 
These maximum energies were not strongly dif- 
ferent in the past during the free expansion stage 
of the remnant evolution since the higher shock 
velocity was almost compensated by the smaller 
shock radius in Eq. (21). 

The similar maximum energies are predicted 
for the core collapse IIP supernovae. They have 
large ejected masses about Mej ^ IOMq and cor- 
respondingly relatively low expansion ve locitie s of 
the order of 3000 km s^^ (see Chevalier (|2005l ) for 
a review). 

The situation is different and the maximum 
particle energies can be higher in the case of the 
core collapse Ib/c and lib supernovae. These su- 
pernovae have high initial velocities about 3 • 10** 



km 



small ejected masses Me 



1 -^ 3 Af, 



o 



and the circumstellar medium corresponding to 
the dense red supergiant stellar wind (lib super- 
novae ) or the interaction zone between the fast 
stellar wind of a Wolf-Rayet progenitor and the 
slow red supergiant wind (Ib/c supernovae) with 
strong magnetic fields (see Chevalier ([2005) for a 
review). In the present paper we have considered 
the acceleration of particles and the generation of 
the MHD turbulence at the parallel shock and our 
theory can not be directly applied to this type of 
supernovae because the magnetic field in the stel- 
lar wind is azimuthal. However, it is possible that 
stellar winds contain significant random magnetic 



fluctuations and we may expect that some part of 
a SNR shock surface may be treated as a paral- 
lel shock. This also will provide the injection of 
thermal particles into the diffusive shock acceler- 
ation process since it is known that the injection 
of thermal ions occurs pre ferenti ally at the paral- 
lel shocks (see Volk et al. (|2003l ) for discussion of 
these topics). 

The Cas A is a lb type supernova. The max- 
imum energy was larger in the past because the 
small radius of the shock in Eq. (21) was com- 
pensated by the higher gas density of the stellar 
wind and the higher shock velocity. Thus the ve- 
locity Ml = 10** km gives the maximum energy 
PmC = 1.2 Z PeV. Therefore such supernovae can 
accelerate cosmic ray protons up to PeV energies. 

Because the maximum momentum (21) de- 
creases at the Sedov stage when the age of the 
remnant increases, the highest energy particles 
leave the remnant. The overall spectrum produced 
by the SNR is form ed in this manner (Ptuskin 
& Zirakashvili (|2005f l) with the energy spectrum 
close to E~'^. At the earlier free expansion stage 
when only a small amount of the supernova ejecta 
energy is transferred to the supernova shock the 
steep high energy t ail in the spectrum is formed 
( Bere zhko & Volk (|2004l ). Ptuskin & Zirakashvili 
(l2005h 'l. This means that the particles from the 
end of the cosmic ray spectrum produced by the 
given supernovae are accelerated at the end of the 
free expansion stage when the shock velocity is of 
the order of the characteristic ejecta velocity. It 
is about 7000-=- 10000 km s"! for la/b/c and lib 
supernovae. 

We should note that the shock velocities in Ta- 
ble 2 are based on radio and X-ray expansion mea- 
surements. If we use the lower shock velocities 
from Volk et al. (|2005h . the corresponding maxi- 
mum energies are a factor of 2 smaller. 

It is clear from Fig. 7 that the amplified mag- 
netic field does not drop downstream of the shock. 
Its level is maintained by the turbulent motions 
produced by the interaction of density inhomo- 
geneities with the shock front. However, the am- 
plitude of these motions slowly decreases down- 
stream of the shock. At some distance magnetic 
and kinetic energies will become equal to each 
other. At larger distances magnetic dissipation 
may occur. These distances are larger than 2-^3^^ 
according to Fig. 7. We conclude that the dissipa- 
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tion length of the magnetic field downstream of 
the shock is not smaller than O.IL in our numeri- 
cal simulation. 

The interaction of the shock front with density 
disturbances results in the shock front deformation 
(see Fig. 6). The thickness of X-ray rims produced 
by the synchrotron cooling of accelerated electrons 
should increase correspondingly up to the values 
about O.OIL according to Fig. 7 (a so-called projec- 
tion effect is not taken into account here). Proba- 
bly this is the reason for the relatively small value 
of the magnetic field 485 /iG in Gas A inferred 
from the width of X-ray filaments in comparison 
with the theoretical value expected at rj^sc > 0.05, 
see Table 2. 

The characteristic time of the random motions 
of the shock front is given by the ratio of the size 
of density inhomogeneities and the shock velocity. 
It is about one year for the size of density inhomo- 
geneities 10^^ cm in young historical SNRs. The 
value L = 0.5R = 1.5 pc was assumed for this 
estimate. 

PDF of the magnetic field downstream of the 
shock that is shown in Fig. 9 is described with a 
good accuracy by the following function: 



Pb{B) 



exp (-V6B/{B^) 



1/2 



(23) 



The exponential tails of the magnetic PDF appear 
to be a fairly universal feature of turbulent ly am - 
plified magnetic fields (Brandenburg et al. ( 19961 ). 
Schekochihin et al. (|2004l )). 

As one can see from Table 1 , the electric poten- 
tial (p has rather large values for the fast shocks 
with velocities u[ > 30000 km s"!. The mean 
electric field is directed opposite to the direction 
of the diffusive electric current. This electric field 
drags the particles which produce the instability 
toward the shock. It may create the steepening of 
the spectrum of these particles. If there is a small 
amount of oppositely charged particles (electrons), 
their spectrum should be somewhat flatter. 

7. Conclusion 

We have investigated the acceleration of par- 
ticles at the fast plane parallel shock. The gen- 
eration of MHD turbulenc e by t he non-resonant 
streaming instability (Bell ( 20041 )) was taken into 



account. We solved this problem using only first 
principles and with a minimum of simplifying as- 
sumptions. We combined the analytical solution 
for the particle acceleration and the numerical 
MHD calculations for the evolution of the MHD 
turbulence. The following results were obtained: 

1) For the relatively fast shocks when the con- 
dition (22) is satisfied, the particles at the high- 
energy end of the spectrum are scattered by small- 
scale random magnetic fields generated by the 
non-resonant streaming instability. Their spec- 
trum has the universal shape given by Eq. (8). 

2) The MHD turbulence is mainly generated by 
the streaming of run-away particles at large dis- 
tance from the shock. The level of the MHD tur- 
bulence is the highest in the shock vicinity. The 
accelerated particles are concentrated in the same 
region. This means that the acceleration of parti- 
cles can be considered in the one-dimensional ap- 
proximation even for a three-dimensional system. 
The characteristic width of the particle distribu- 
tion in our simulations is not larger than 0.1 i (see 
Fig. 3 and 10th line of Table 1). 

3) It is important, that the non-resonant in- 
stability produces strong density fluctuations up- 
stream of the shock (see Sect. 4 and Paper I). These 
fluctuations produce a strong deformation of the 
shock front and fast vortex motions downstream 
of the shock. That is why the magnetic amplifica- 
tion in the shock transition region is not reduced to 
a simple compression of the magnetic field in the 
direction perpendicular to the shock front. The 
magnetic field is also stretched by the fiow motions 
in the direction perpendicular to the shock front. 
As a result, the magnetic field component which is 
perpendicular to the shock front is a factor of 1.4 
larger than the parallel components downstream 
of the shock. This naturally explains the prefer- 
able radial orientation of magnetic fields in young 
SNRs. 

The characteristic time 1 year of the shock de- 
formation and of the corresponding MHD fluctu- 
ations found here (see Sect. 6) is of particular in- 
terest in the light of the last results on variabil- 
ity of X-ray emission observed in RXJ1713 SNR 
(Uchiyama et al. [200^). 

4) The dissipation of the magnetic held down- 
stream of the shock is relatively slow in our simula- 
tion. If so, the origin of X-ray fllaments, observed 
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in young SNRs, is related to the fast synchrotron 
coohng of accelerated electrons but not to the de- 
cay of the MHD turbulence. 

5) The values of the calculated amplified mag- 
netic field are similar to those, observed in his- 
torical SNRs, if the energy flux of the run-away 
particles is not low: r/esc > 0.05 (see Table 2). 

6) The magnetic field growth is only linear in 
time for the fast shocks with the normalized ve- 
locities higher than about ten thousand km s~^. 
This reduces the maximum energy of accelerated 
particles compared to the case of the exponential 
growth. For these shocks the energy of the mag- 
netic field amplified upstream is a small fraction 
s/Va/c of the energy density of the highest energy 
particles (see Eq.(20)). The magnetic field ampli- 
fication is relatively weak for slow shocks with the 
normalized velocities 1230 km s~^ and smaller (see 
Eq.(18)). 

7) We calculated numerically the maximum en- 
ergy of accelerated particles (see Fig. 2). This 
energy may be described by the analytical for- 
mula (19). The maximum energies of particles 
are higher than the energies obtained in the Bohm 
limit in the background magnetic field (see the 7th 
line of the Table 1) but lower than the energies ob- 
tained using the Bohm limit in the amplified field. 
The significantly long time of the random mag- 
netic field growth is the main factor that limits 
the maximum energy of accelerated particles. We 
calculated the maximum energies for four histori- 
cal SNRs (Table 2). 

8) Using the last result we found that the max- 
imum energy of cosmic ray protons accelerated by 
la and IIP supernovae is about 100 -j- 300 TeV. 
Only Ib/c and presumably lib supernovae may 
accelerate protons up to PeV energies. Since it 
is expected that the explosion rate is the highest 
for IIP supernovae, we should observe a change of 
the slope of the galactic cosmic ray spectrum at 
the energies of the order of 100 TeV. 

9) The MHD turbulence generated by the non- 
resonant streaming instability has a non-zero mag- 
netic helicity. The helical Lorentz force produces 
corresponding plasma motions and the mean elec- 
tric field Eq that is in the opposite direction to the 
electric current of energetic particles. This electric 
field modifies the cosmic ray transport equation 
(see Paper I) . This effect is significant for very fast 



shocks with velocities larger than 30 thousands km 
s~^. The presence of this field should result in 
the steepening of the spectrum of particles which 
produce the non-resonant instability (presumably 
nucleons) and the flattening of the spectrum of op- 
positely charged particles (presumably electrons). 
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of valuable suggestions. We are grateful to Heinz 
Volk for many fruitful discussions of acceleration 
by astrophysical shocks. VNZ and VSP acknowl- 
edge the hospitality of the Max-Planck- Institut fiir 
Kernphysik, where this work was mainly carried 
out. The work was also supported by the RFBR 
grant in Troitsk. 
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A. The dependence of scattering on pitch angle. 

If the random field is not isotropic the same is true for the particle scattering. Let us assume that the 
random magnetic field is isotropic in the plane perpendicular to the mean magnetic field Bq and that the 
cosmic ray distribution /o depends only on pitch angle 0. Now the scattering operator described by the 
tensor Uij (cf. Paper I) can be written as 

where /i — cob 9. The scattering frequency vo{p-) can be expressed in terms of the spectrum of y-component 
of the random magnetic field (mean field is in z direction): 

vo(fi) = 27r-^ / cPkByJ\i)6{k^n + sin6'). (A2) 
P J 

The parallel diffusion coefficient D\\ can be written now as (see e.g. Berezinskii et al. ( 1990l) ) 

^11 - T / (A3) 

2 J vq{h) 
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